OPEN 3 ACCESS Freely available online 



•0-PLOS I ONE 



Characterization of Human Papillomavirus Type 154 and 
Tissue Tropism of Gammapapillomaviruses 

Agustm Enrique Ure, Ola Forslund* 

Department of Laboratory Medicine, Section of Medical Microbiology, Lund University, Malmo, Sweden 



Abstract 

The novel human papillomavirus type 1 54 (HPV1 54) was characterized from a wart on the crena ani of a three-year-old boy. 
It was previously designated as the putative HPV type FADI3 by sequencing of a subgenomic FAP amplicon. We obtained 
the complete genome by combined methods including rolling circle amplification (RCA), genome walking through an 
adapted method for detection of integrated papillomavirus sequences by ligation-mediated PCR (DIPS-PCR), long-range 
PCR, and finally by cloning of four overlapping amplicons. Phylogenetically, the HPV154 genome clustered together with 
members of the proposed species Gammapapillomavirus 11, and demonstrated the highest identity in LI to HPV136 
(68.6%). The HPV154 was detected in 3% (2/62) of forehead skin swabs from healthy children. In addition, the different 
detection sites of 62 gammapapillomaviruses were summarized in order to analyze their tissue tropism. Several of these 
HPV types have been detected from multiple sources such as skin, oral, nasal, and genital sites, suggesting that the 
gammapapillomaviruses are generalists with a broader tissue tropism than previously appreciated. The study expands 
current knowledge concerning genetic diversity and tropism among HPV types in the rapidly growing gammapapillo- 
mavirus genus. 
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Introduction 

The papillomaviruses (FVs) are small viruses with icosahedral 
symmetry and a circular, double stranded genome [1]. These 
viruses are widely distributed across vertebrates, and among 
humans alone, more than 170 unique complete genomes of types 
have so far been sequenced [2]. The papillomaviruses are 
epitheliotropic and produce hyperproliferation of squamous cells 
known as papillomas. Human papillomaviruses (HPV) are divided 
into high-risk types, which are etiologically associated with cancer 
of the cervix uteri; and low-risk types, which produce benign warts 
[3,4] . Regarding the taxonomy of papillomaviruses, the Interna- 
tional Committee for the Taxonomy of Viruses (ICTV) recom- 
mends inclusion of both genotypic and phenotypic information in 
the definition of viral genera and species [5]. Nevertheless, 
papillomaviruses have been an exception to the classical rules, and 
a system based mainly on sequence identity was adopted (reviewed 
in de Villiers, 2013). Accordingly, a 70% sequence identity 
amongst LI ORFs is used as cut-off guide to classify viruses as 
belonging to the same species [6]. Even so, there is a gray zone 
(67.5%-70.5% identity) where the distribution curves for inter- 
and intraspecies identities overlap, and consequently the classifi- 
cation must be curated [7] . The papillomaviruses have in the past 
been classified as mucosal or cutaneous according to their tropism. 
The human cutaneotropic papillomaviruses are represented 
mainly by the genera Beta- and Gammapapillomavirus ((3-PV and y- 
PV). Their prevalence and natural history, for example acquisition 
soon after birth, reflects a commensalic interaction with the 



immunocompetent host [8-10]. Historically, several cutaneous 
papillomaviruses have been isolated from patients with the genetic 
disease epidermodysplasia verruciformis [11,12], and from immu- 
nosuppressed patients [13], probably due to higher viral loads 
among these patient groups [14]. As a consequence of improved 
methods, several PVs have been characterized, and the number of 
HPV types of the y-PV genus has been expanded from 1 6 HPV 
types in 2010 [7] to 62 representative HPV types in 2013 
(retrieved from GenBank, September 2013). A review of the 
isolation sources of these HPV types could lead to increased 
knowledge of their tropism. Here, we categorized the different 
isolation sites of the y-PV as genital, oral, nasal or cutaneous 
(including skin lesions and healthy skin). The aim of the study was 
to obtain the complete sequence of a novel papillomavirus, and to 
present a summary of the isolation sources of the rapidly growing 
genus Gammapapillomavirus. 

Results 

HPV 1 54 was isolated from a wart on the crena ani of a three- 
year-old boy. The index sample was negative for HPV by PCR 
using MGP primers [15] and a Luminex system [16], but positive 
for HPV by FAP-PCR [10]. The FAP amplicon was sequenced 
and showed the highest identity to the FADI3 fragment (99.3%, 
acc. no. FJ480954). The FADI3 was originally amplified from a 
forehead swab of a six-year-old girl [17]. 

The FAP amplicon represented a putative novel type, as the 
partial LI ORF was below 90% of sequence identity to any other 
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known PV type [6]. With or without pre-amplification by rolling 
circle amplification (RCA), we railed to amplify the complete 
genome by long-range PCR or with PCRs designed to cover half 
the genome by a combination of LI -specific primers with 
degenerated El primers (data not shown). As a consequence, we 
adapted the method for detection of integrated papillomavirus 
sequences by ligation-mediated PCR (DIPS-PCR) [18] to obtain 
sequence information outside the FAP amplicon region of the LI 
ORF. In order to test the performance of the DIPS-PCR method, 
we verified the sequence around the integration site of HPV16 into 
the genome of SiHa cells (data not shown). Several DIPS-PCRs 
and PCR with degenerated HPV primers were used to obtain 
sequence data from L2 to almost the end of LI (Figure 1). 
Proximal to the ends of that region, new primers were designed 
that were combined with degenerated primers based on related y- 
PVs, and two amplicons of ~2000 bp were cloned and sequenced. 
In order to obtain the complete genome, four overlapping 
amplicons were obtained by PCR with specific primers 
(Figure 1). The viral load of HPV154 was 68 genomes per human 
cell (Table 1). 

The four clones were submitted to the International Reference 
Center for Human Papillomaviruses at the German Cancer 
Research Center, Heidelberg, Germany, and the compiled 
sequence was verified and officially designated as HPV 154 
(GenBankJN211193). 

The genome of HPV154 comprised 7286 bp with a GC content 
of 37.9%. The most closely related type was HPV136, isolated 
from the oral cavity [19], with 68.6% identity at the LI ORF, 
whereas the uncloned HPV isolate KN3, obtained by high 
throughput sequencing from healthy skin [20], showed 71.8% 
identity at LI. 

HPV 154 demonstrated the classical genomic organization of 
PVs, with seven ORFs identified (Figure 1). The putative E6 
protein had two zinc-finger domains (CX 2 CX 2 9/3oCX 2 C) that 
were separated by 36 amino acids, which are conserved among 
PVs [21]. Similarly, there was one zinc -finger domain in the 
inferred E7 protein, as well as the tumor suppressor (pRB) binding 
domain (LXCXE) [22,23]. In the putative El protein, the 
superfamily 3 ATP-dependent helicase domain was identified 
[24,25]. The theoretical E2 protein had the typical C-terminal 
DNA-binding domain and the N-terminal trans-activation domain 
[26,27]. The E4 ORF showed a start codon; nevertheless it was 
ignored as we identified the characteristic donor (AAG/ 
GUASNR) and acceptor (GUYACYAG/YU) RNA splicing sites 
[28], and the resulting putative E1 A E4 fusion protein with six 
amino acids from the El N-terminal end. 

The early polyadenylation site (AATAAA) for processing of 
early mRNAs, was located at the 5' end of the L2 ORF, while the 
late polyadenylation site was downstream of LI within the 
upstream regulatory region (URR). The URR was relatively 
short, being 517 bp, and contained six putative E2 binding sites 
(E2BS, ACCN 6 GGT). Near the E2BS proximal to E6, several El 
binding sites (AACAAT) or related AT-rich sequences were 
identified and probably represent the origin of replication [29,30]. 
A TATA box was found (pos. 6973-6977) and surrounded by two 
E2BS in close proximity (2 and 14 bp). 

As the FADI3 amplicon and HPV154 were obtained from 
samples from children, we tested whether this PV type might have 
an increased presence in this group. Among the cutaneous swab 
samples, HPV154 was detected in 3% (2/62) by real-time PCR. 
One of the samples showed a viral load of 0.093 genomes and the 
other 1.3 genomes per human cell (Table 1). 

The phylogenetic relationships of HPV 154 were inferred based 
on the complete genomes of 91 HPV-sequences (Figure 2). 



HPV 154 was positioned on a divergent group, along with the 
members of the currently proposed y-PVll species (pending 
ICTV approval) (Figure 2). HPV154 is below the 70% identity 
limit suggested for species definition compared to its closest 
relative type, HPV136 (68.6%). Moreover, the lowest values of 
identity for HPV 154 within y-PVll (Figure 2) are with HPV 140 
(65.37%) and HPV169 (65.86%). 

In order to increase knowledge of the biological niches of 
gammapapillomaviruses, we enumerated the isolation sources 
(complete HPV-genomes) and detection sites of all HPV types 
from this genus (Table 2). The y-PVl species contains members 
isolated from a wide variety of sites, including healthy tissue such 
as the skin, nasal cavity and male genitalia. However, they are also 
found in non-healthy tissue such as common warts, pigmented 
verrucas, squamous cell carcinoma (SCC) and actinic keratosis 
(AK). In a similar fashion, HPV types of the y-PV3, y-PV7, y- 
PV8, y-PV9, y-PVIO, y-PVll and y-PV15 species have been 
found in several different sites, including both mucosal and 
cutaneous sources (Table 2). The y-PV6 members (HPV101, 103 
and 108) are the only y-PV restricted to genital and oral mucosal 
sites (six isolates). On the other hand, the y-PV 12 members have 
been isolated only from warts and healthy skin (eight isolates). The 
remaining species, which include y-PV2, y-PV4, y-PV5, y-PV 13, 
y-PV 14, y-PV 16, y-PV 17 and the six putative novel species 
(Figure 2), have members with less than three findings, but include 
both mucosal and cutaneous sites (Table 2). 

Discussion 

From a wart on the crena ani of a three-year-old boy we 
characterized the novel HPV 154, which was added to the rapidly 
growing genus Gammapapillomavirus. The genome of the HPV 154 
demonstrated the typical genomic organization of papillomavirus- 
es and was lacking the E5 ORF, as with y- fi- and u-HPVs [31]. 
The prevalence of HPV 154 was 3% in skin swabs from healthy 
children, which is in a similar range as other HPV types of the beta 
and gammapapillomavirus genera detected among adults [32-38]. 
A limitation of the study was that only samples from children were 
used in the screening of HPV154. Another limitation was that 
HPV 154 was characterized from a swab sample, so we cannot 
assert that the virus caused the lesion. Our initial approach of 
obtaining the complete HPV genome by long-range PCR failed, 
which may be due to the low concentration of HPV154 (Table 1) 
and/or fragmentation of the genome. Instead we used PCRs with 
degenerated HPV primers and an adapted version of the DIPS- 
PCR method [18]. It was notable that several sequence reads of 
the restriction enzyme site for the adapter ligation of DIPS-PCRs 
could not be verified in the final sequence, which might indicate 
star activity of the restriction enzyme prior to ligation of adapters. 
Thus, compiled sequences from DIPS-PCR should be verified by 
independent PCR. Nevertheless, we showed that the DIPS-PCR 
can expand the length of sequence information, which in turns 
allows to design HPV primers at positions that reduce the size of 
the targeted long-range amplicons, which would probably increase 
amplification efficiency. However, our approach is laborious and 
time-consuming compared to modern methods such as high- 
throughput sequencing, and would be a feasible alternative only in 
cases where the latter methods are not available or other classical 
methods such as RCA and long-range PCR have failed. 

HPV 154 demonstrated the typical genomic organization of 
papillomaviruses, but lacked the E5 ORF as in all known members 
of the |3- y- and |X-PV genera [31]. The prevalence of HPV154 
(5%) in skin swabs from healthy children was in a similar low range 
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URR clone 



E clone 




L clone 



2 clone 



Figure 1. Methodological strategy and genomic organization for HPV154. The ORFs are indicated with light blue arrows. The primers used 
in the DIPS-PCR are represented with half arrows. The inner circle shows the strategies employed, and the outer circle shows the final clones. Putative 
binding sites for viral proteins and cellular factors are shown as follows: E2 binding site, E2BS (T); El binding site, E1BS (A); TATA-box (■); 
Polyadenylation signal (•). 
doi:1 0.1 371 /journal.pone.0089342.g001 



as for other HPV types of the (3- and y-PV genera detected among 
adults [32-37]. 

Concerning the presented phylogenetic tree, HPV 154 was 
placed with high support along with members of the currently 
proposed y-PV 11 species [2]. Nevertheless, the percentages of LI 
ORF sequence identity among different members of this group 
were closer to the interspecies center of the distribution (64.5%) 
than to the intraspecies center (72%) as shown in the histograms of 



LI sequence identity from 189 PVs [7]. Even using 68.5% of 
identity as a more stringent/ divergent limit for inter-species 
distance, three different species could be delimited, as shown with 
a gray dotted key in Figure 2. 

Accordingly, a more stringent/divergent limit of 68.5% identity 
allowed us to split the proposed y-PVl 1 into three putative novel 
species. However, the proposed y-PVl 1 [2] were based on HPV 
types with tropism for the oral cavity, but additional HPV types 
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Table 1. Characteristics of HPV154 positive children and viral loads in swab samples. 





Sample ID 


Age 


Gender 


Swab sampling site 


Sample type 


Viral genomes/nL (CV%)* 


Viral genomes/human cell (CV%)* 


1* 


3 


Male 


Wart, crena ani 


Extracted DNA 1 ' 


48 (6.3) 


68 (7.7) ¥ 


2 


3 


Male 


Wart, crena ani 


Saline 


41 (20) 


39 (95) 


3 s 


1 


Female 


Forehead 


Saline 


0.49 (73) 


0.093 (59) 


4 5 


4 


Female 


Forehead 


Saline 


0.68 (69) 11 


2.0 (23)'' 



*lndex sample used for characterization of HPV154. 

1 DNA was extracted from sample no. 2. 

"'"Viral load measurements were performed in triplicate. 

¥ Number of human cells obtained in duplicate due to one unsuccessful measurement. 
§ HPV154 detected in 2 forehead swab samples from 62 children. 

^'Number of HPV154 copies obtained in duplicate due to one unsuccessful measurement. 
doi:1 0.1 371 /joumal.pone.0089342.t001 



within this species had been isolated from healthy skin and warts as 
in the case of HPV154 (Table 2). According to ICTV "A species is 
a monophyletic group of viruses whose properties can be 
distinguished from those of other species by multiple criteria" 
(http://www.ictvonlfne.org/ codeOfVirusClassification.asp), but 
among a majority of the y-PV species it is difficult to find 
distinctive characteristics besides sequence identity. Although 
species definition for y-PV is currently based mainly on sequence 
identity [2,6,7], the criteria have varied over time, and hence have 
not been applied homogeneously to all members of the genus. For 
example, HPV48 (y-PV2) shares 68.01% identity with HPV50 (y- 
PV3) and 69.92% with HPV131 (proposed y-PV14). A system 
with well supported phylogeny would help to define species, and it 
may be safer to leave new members unclassified at the species level 
until more isolates are sequenced. 

With regard to the isolation sites enumerated in Table 2, it 
seems that y-PVs can be found in a broader variety of sites than 
previously appreciated. Diverse sites are found for members of the 
same species or even type. Even though very few studies have 
searched for papillomavirus in the oral or nasal cavities [19,39-41], 
members belonging to nearly every species have been identified in 
those studies. Altogether, this supports the idea that y-PVs are 
generalist with various forms of tropism. The only exceptions seem 
to be the members of the y-PV6 species, HPV101, HPV103 and 
HPV108, as they appear to have only mucosal tropism. This group 
also distinctively lack the E6 ORF and have been associated with 
disease [42,43]. 

However, our approach to summarize isolation sites of 
gammapapillomaviruses described in reports provides only a 
suggestive mode to study the tropism of these viruses. In order to 
perform a rigorous evaluation of the tropism for each HPV type, a 
random population study with samples collected from several sites 
is needed. 

In this study we successfully characterized the novel HPV 154, 
thereby expanding knowledge about the diversity and tropism of 
HPV types in the rapidly growing y-PV genus. In addition, we 
suggest that gammapapillomaviruses are generalist with broad 
tissue tropism. We expect that improved amplification methods 
will promote the discovery of additional HPV types of the y-PV 
genus, and will shed light onto its apparent wider tropism 
compared to other papillomavirus genera. 

Materials and Methods 

Ethics Statement 

Ethical approval for this study was granted by the Ethical 
Committee of Lund University (LU 106-01). The samples were 



obtained with written informed consent from the parents of the 
minors. 

Sample Processing: DNA Extraction 

A swab suspension of 200 |tl was processed using the automated 
MagNA Pure LC with the Total Nucleic acid kit (Roche), and 
eluted in 100 ul. 

Standard HPV Analysis 

Sample adequacy was assessed by testing 5 |xL of the sample for 
the human fi-globin gene with a real-time PCR [44]. For 
identification of 22 genital HPV types (6, 11, 16, 18, 31, 33, 35, 
39, 42, 45, 51, 52, 56, 58, 59, 66, 68, 70, 73, 82) simultaneously 
5 |lL of extracted material was added to a total volume of 25 |iL 
for MGP-PCR and subsequent Luminex analysis [15,16]. Five 
microliters was used for FAP-PCR targeting the LI ORF [10]. 

RCA 

The eluted total DNA was subjected to rolling-circle amplifi- 
cation (RCA) using illustra TempliPhi 100 Amplification Kit (GE 
Healthcare), basically following the manufacturer's instructions, 
but in a slightly modified form. As indicated elsewhere [45], the 
final concentration of each dNTP was increased to 450 uM and 
the reaction was incubated overnight. 

Genome Walking with DIPS-PCR 

The RCA product (~0.5 |ig) was digested in order to be used as 
input for DIPS-PCR [18] with either TaqI, Sau3AI (Bspl43I), 
FatI, Xbal or Hindlll enzymes (10 units, 20 J_ll final volume) 
following the manufacturer's instructions (Fermentas). This 
method employs short adapters specific to each enzyme. These 
short oligos were synthesized to be compatible with enzymes TaqI 
(cgcaacgtgtaagtctg), Sau3AI (gatccaacgtgtaagtctg), FatI (catg- 
caacgtgtaagtctg), Xbal (ctagcaacgtgtaagtctg), and Hindlll (agct- 
caacgtgtaagtctg); and all of them were modified to contain a 
phosphate group at the 5' end and an amino group (Amino C3) at 
the 3' end (MWG, Germany). The last three enzymes were added 
in this study in order to increase the likelihood of obtaining longer 
fragments and to speed up the overall process. Each of the short 
adapters was combined with a unique long adapter (gggccatcagt- 
cagcagtcgtagccggatccagacttacacgttg) in equimolar quantities, to a 
final concentration of 25 (J.M. The combined oligos were 
incubated in a thermocycler for 4 min at 94°C, and then 
incubated for 5 s at 95°C. This process was repeated 18 times, 
reducing the temperature by 5°C in each step. The formed 
adapters were then aliquoted and stored at — 20°C. The RCA 
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Dyoxipapillomavirus 
Pipapillomavirus 




Figure 2. Phylogenetic tree for HPV154 and related papillomaviruses. Ninety one complete genomes were analyzed, including all 
gammapapillomaviruses as well as related genera but closer than betapapillomaviruses. HPV1 54 is indicated with an arrow. The tree was obtained by 
the maximum likelihood approach using RAxML software and rooted with the betapapillomavirus HPV5 and HPV9. Bootstrap support values are 
indicated in each branch as percentages. Genera are indicated at the rightmost side, while species of y-PVs are to the right of brackets. The dotted 
key to the left of 7-PVH* shows species suggested in this work. The bar indicates the number of substitutions per site. * Proposed category and 
pending approval by ICTV. fPutative novel species, Uncloned HPV genome (not official type). The following species are represented by the 
types included in the tree: AsPVl, Apodemus sylvaticus papillomavirus 1; BPV, Bos taurus papillomavirus spp.; ChPV1, Capra hircus papillomavirus I; 
CPV, Canis familiaris papillomavirus spp.; DdPV1, Delphinus delphis papillomavirus i; HPV, Gammapapillomavirus spp.; McPV2, Mastomys coucha 
Papillomavirus 2; PphPV, Phocoena phocoena papillomavirus spp.; PsPV1, Phocoena spinipinnis Papillomavirus I; TtPV, Tursiops truncatus Papillomavirus 
spp. 

doi:1 0.1 371 /journal.pone.0089342.g002 
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Table 2. Sites of detection of gammapapillomaviruses. 



HPV type 


Species 


Site 


References 


FA-fragment 




(% identity to type) 


HPV4 


y-PVl 


Common wart 


[64] 










Male genitalia 


[57] 










Nasal cavity 


[39] 










Skin 


[62] 






HPV65 


y-PVI 


Pigmented verruca 


[69] 










SCC/AK 


[38] 






HPV95 


y-PVI 


Skin wart 


[71] 






HPV48 


Y-PV2 


sec 


[72] 






HPV50 


Y-PV3 


Skin macules EV 


[56] 








" 


Nasal cavity 


[39] 










Skin 


[9,61,62] 








" 


SCC/AK and keratoacanthoma 


[38] 






HPV60 


Y-PV4 


Planta cyst 


[63] 






HPV88 


Y-PV5 


Skin wart 


[36] 






HPV101 


Y-PV6 


Cervicovaginal 


[42] 










Male genitalia 


[57] 






HPV103 


Y-PV6 


Cervicovaginal 


[42] 






" 


" 


Male genitalia 


[57] 










Nasal cavity 


[39] 






HPV108 


Y-PV6 


Cervicovaginal 


[43] 






HPV109 


Y-PV7 


sec 


[37] 










Seborrhoeic keratosis 


[37] 










Cervicovaginal (CINI) 


[37] 










sec 


[65] 


FA137 (92.7) 




HPV123 


Y-PV7 


Oral cavity 


[19] 










Skin SCC/AK 


[38] 










Skin Seborheic keratosis 


[65] 


FA136 (94.7) 




HPV134 


Y-PV7 


Oral cavity 


[19] 










Male genitalia 


[57] 










Nasal cavity 


[39] 










AK 


[10] 


FA1.1 (94.6) 








Skin 


[10] 


FA 1.1 (94.6) 




HPV138 


Y-PV7 


Oral cavity 


[19] 








" 


Nasal cavity 


[39] 










Skin 


[8] 


FA20 (91.5) 




HPV139 


Y-PV7 


Oral cavity 


u. d. 










Cutaneous adnexal tumour 


[66] 


FA105 (92.2) 




HPV149 


Y-PV7 


Skin wart 


[13] 






HPV155 


Y-PV7 


Actinic keratosis 


[38] 










Nasal cavity 


[39] 


FA166 (91.4) 




HPV170 


Y-PV7 


Skin 


[59] 






HPV112 


Y-PV8 


Condyloma acuminata 


[37] 






HPV119 


Y-PV8 


Oral cavity 


[19] 










Male genitalia 


[57] 










Keratoacanthoma 


[38] 










Skin 


[67] 


FA155 (90.8) 




HPV147 


Y-PV8 


Oral cavity 


u. d. 










Nasal cavity 


[39] 
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Table 2. Cont. 





HPV type 


Species 


Site 


References 


FA-fragment 


(% identity to type) 






Skin 


[10] 


FA13 (94.7) 


HPV164 


Y-PV8 


Skin 


[59] 




HPV116 


Y-PV9 


Rectal Swab 


[68] 




HPV129 


Y-PV9 


Skin wart 


[13] 




FD2* 




Skin 


u. d. 




HPV121 


Y-PVIO 


Oral cavity 


[19] 








Male genitalia 


[57] 




HPV130 


Y-PVIO 


Skin wart 


[13] 








Male genitalia 


[57] 




" 




Nasal cavity 


[39] 








Actinic keratosis 


u. d. 


FA147 (95.1) 


HPV133 


y-PVIO 


Skin wart 


[13] 






" 


Male genitalia 


[57] 








Nasal cavity 


[39] 








Skin 


[61] 


FA44 (92.4) 


HPV142 


Y-PVIO 


Oral cavity 


u. d. 








Nasal cavity 


[39] 








Skin 


[9] 


FA67 (92.2) 


FA69* 


Y-PVIO 


Condyloma acuminata 


[60] 








Skin 


[9] 


FA69 (100) 


HPV126 


Y-PV11* 


Skin wart 


[70] 




HPV136 


Y-PV11* 


Oral cavity 


[19] 








Skin 


[10] 


FA8 (91.8) 


HPV140 


Y-PV11* 


Oral cavity 


u. d. 




HPV141 


Y-PV11* 


Oral cavity 


u. d. 




HPV154 


Y-PV11* 


Skin wart 


Present study 








Skin 


[17] 


FADI3 (91.3) 


KN3* 




Skin 


u. d. 




HPV169 


Y-PV11* 


Skin 


[59] 




HPV132 


Y-PV12* 


Skin wart 


[13] 








Skin 


[61] 


FA78 (91.3) 


HPV148 


Y-PV12* 


Skin wart 


[13] 




HPV165 


Y-PV12* 


Skin 


[59] 




CG2* 


Y-PV12* 


Skin 


[20] 








Skin 


[66] 


FA106 (96.2) 


CG3 1 ' 


Y-PV12* 


Skin 


[20] 




KC5* 




Skin 


[59] 




HPV128 


Y-PV13* 


Skin wart 


[13] 




Mrv I dd 


D\/1 3* 

y-r V 1 3 


Condyloma acuminata 


L44J 




HPV131 


Y-PV14* 


Skin wart 


[13] 




HPV135 


Y-PV15" 


Oral cavity 


[19] 








Male genitalia 


[57] 








Nasal cavity 


[39] 








Skin 


[8] 


FA28 (97.5) 


HPV146 


Y-PV15* 


Oral cavity 


u. d. 








Nasal cavity 


[39] 








Skin 


[9] 


FA35 (96.8) 


" CGI* 




Skin 


[20] 
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Table 2. Cont. 





HPV type 


Species 


Site 


References 


FA-fragment 


(% identity to type) 


SD2* 




Nasopharyngeal and oropharyngeal 


[73] 




HPV137 


Y-PV16* 


Oral cavity 


[19] 




HPV144 


Y-PV17* 


Oral cavity 


[19] 








Male genitalia 


[57] 




HPV156 


v-PV^ 
i > v 


Skin 


[58] 




Fi864* 


y-pyt 


Nasal swab 


[40] 




HPV161 


y-PV 1 " 


Skin 


[59] 




HPV162 




Skin 


[59] 




HPV166 




Skin 


[59] 




KN1* 


y-pyt 


Skin 


u. d. 




FD1* 




Skin 


u. d. 




FS1 * 




Skin 


u. d. 




SE87* 




Condyloma acuminata 


[60] 




KN2* 


y-PVt 


Skin 


u. d. 





*Proposed species pending on ICTV approval. 
1 Putative novel species. 

"'"Complete genomes, not yet designated HPV numbers. SCC, Squamous Cell Carcinoma; AK, Actinic keratosis; EV, Epidermodysplasia verruciformis; u. d., unpublished 
data. The first row for every HPV type represents the reference clone (complete genome). 
doi:1 0.1 371 /joumal.pone.0089342.t002 



product previously digested was combined with the appropriate 
(enzyme specific) double stranded adapter (50 pmol) and ligated 
by adding ATP (0.5 mM), DTT (10 mM) and 1.75U of ligase 
(Fermentas) to a final volume of 2 7 jiL The reaction was incubated 
at 16° overnight, and later diluted to 40 (J.1 with water and was 
used as template for all subsequent PCR steps. Two microliters of 
the ligation products were subjected to an initial linear PCR 
amplification using a viral specific primer (0.2 |lM) in lx PCR 
buffer (Roche), 1.8 mM MgCl 2 , 0.2 mM each dNTP, 1.25 U of 
Taq polymerase (AmpliTaq Gold, Roche). The cycling conditions 
were 94°C 5 min; 45 cycles of 94°C 15 s, 50°C 20 s, 72° 3 min; 
72° 5 min. A second PCR (exponential) using 2 of the linear 
PCR as a template was made using a nested viral-specific primer 
(primer sequence available upon request) combined with the 
adapter-specific "API" primer (ggccatcagtcagcagtcgtag) in the 
same conditions as the linear PCR but with 0.5 |J,M of both 
primers. The cycling conditions were 94°C 5 min; 35 cycles of 
94°C 15 s, 55°C 20 s, 72°C 3 min; 72°C 5 min. The PCR 
products were evaluated by electrophoresis, and the longest 
product (from the different restriction enzymes) was selected, 
cloned using a TOPO-TA cloning kit (Invitrogen) according to the 
manufacturer's instructions, and sequenced at MWG (Germany). 
The novel sequence was used to design primers with Oligo 7 
software (Molecular Biology Insights, USA) in order to carry out 
successive steps of linear/ exponential PCR. 



and a final extension of 72° 7 min. The PCR product was purified 
with an Illustra Microspin S-300 HR column (GE Healthcare) and 
sequenced with a nested specific primer. 

A second set of degenerated primers was designed for the E 1 
ORF. The primers 'El gamma Fnew' (gacagtggdatwkddgaagat- 
gaa) and 'El gamma Rnew' (ttcatcttcddmwatdccactgtc), were 
combined with 'FADI3 -2 nes R' and 'FADI3 +2 nes F', 
respectively (0.3 uM final concentration). One microliter of 
RCA-pre-amplified template was subjected to Long Range PCR 
(Expand Long Template PCR System, Roche) amplification in lx 
PCR buffer 2 (Roche), 0.5 mM each dNTP, 0.7 U of Enzyme mix 
(Expand Long Template PCR System, Roche). The cycling 
conditions were 94°C 2 min; 10 cycles of 94°C 15 s, 50°C 1 min, 
68°C 4 min and 30 cycles of 94°C 15 s, 59°C 30 s, 68°C 4 min. 
The product was electrophoresed and the bands excised from the 
agarose gel, purified with a QJAquick Gel Extraction kit (Qiagen) 
and cloned using a TOPO-TA cloning kit (Invitrogen), according 
to the manufacturer's instructions. 

The remainder of the genome was amplified using an Expand 
High Fidelity PCR system (Roche, Mannheim, Germany) with 
2.5 ul of the RCA amplified DNA (diluted 1:100) as a template; 
0.3 [LM of each primer, 3.5 mM MgCl 2 and 1.24 U of enzyme in 
a final volume of 50 joJ. The cycling was 94°C 2 min; 40 cycles 
94°C 15 s, 55°C 20 s, 72°C 2 min increasing 5 s/cycle; 72°C 
10 min. The product was purified and cloned as indicated before. 



PCR with Degenerated Primers and Overlapping 
Amplicons 

Primers were designed in the L2 ORF from a multiple 
alignment of HPV sequences (HPV48, 50, 60, 65, 88, 95, 112 
and 116) related to the known FAP region of HPV154. The 
degenerated primer "L2 gamma F" (tttgratwtgaaaatcccgccttt) was 
used in conjunction with the HPV154 specific primer "PV77 
FADI3 R" (gacctgtgctaccgactccaag). The cycling conditions were 
94°C for 5 min; 40 cycles of 94°C 15 s, 52°C 20 s, 72°C 1 min 



Screening of Samples from Children 

Forehead samples from children from three age groups, one 
month (23), one year (19), and four years old (20), were obtained 
from a previous study [9] and used as templates for the screening. 
The samples were swab suspensions in 0.9% NaCl, and were used 
directly as PCR templates. Primers intended for real-time PCR 
were designed using Oligo 7 (Molecular Biology Insights, USA) 
and tested for possible cross-reactivity with other types using 
Prirner-BLAST [46]. The resulting primers were 'HPV154 5511 
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F' (accgtggtggtcctcttgg) and 'HPV154 5647 R' (tgggtcaaag- 
gaaatgttttgg). The real-time PCR was carried out in an ABI 
7500 Real-Time PCR Instrument (Applied Biosystems). Each test 
tube contained Power SYBR Green PCR Master Mix (Applied 
Biosystems), 0.3 U.M of each primer, and 2.5 |J.l of template in a 
final volume of 25 pi. Pipetting was automated using a Qiagility 
device (Qiagen). The cycling conditions were 95°C 10 min and 45 
cycles 95°C 15 s, 60°C 60 s, followed by melting curve analysis. 
For data analysis the ABI 7500 software v2.0.6 was used and the 
threshold for positivity was automatically calculated. The speci- 
ficity of the HPV154 PCR was also analysed by gel-electrophoresis 
for identification of the expected 137 bp amplicon. Plasmid DNA 
concentration of HPV154 clone L (Figure 1) was quantified using 
a spectrophotometer (NanoDrop ND-1000, Nanodrop Technol- 
ogies, Oxfordshire, UK). The sensitivity of the method was 
evaluated using serial dilutions of the clone L (Figure 1), and 10 
copies were detected in a background of 1 ng of human DNA 
(Sigma-Aldrich, art. D 701 1). The positive HPV154 control had a 
Tm of 75.4C (CV: 0.2%, based on four measurements). 

Viral Load of HPV154 Positive Samples 

The number of viral genomes per cell was quantified by 
carrying out two separate quantitative real-time PCR assays to 
amplify a part of the HPV154 L2 gene and the human P-globin 
gene. For the quantification of HPV154 the primers and PCR 
conditions were identical to that used for screening of samples 
from children as described above. Quantification was extrapolated 
from a linear regression standard curve obtained from serial 
dilutions of 100,000 to 100 copies per PCR of HPV154 plasmid 
DNA (clone L, Figure 1) in a background of 10 ng human 
placenta DNA (Sigma-Aldrich, art. D 7011). The standard curve 
had a slope of -3.4, y intercept of 37.3 and r 2 of 0.99. The PCR 
efficiency calculated from slope was 97.5%. Similarly, in order to 
calculate the number of cells analyzed per sample, the P-globin 
gene was amplified with PC03 and PC04 primers in a 25 |xl PCR 
reaction containing 2.5 uL template [44]. The standard curve was 
obtained from serial dilutions of 50,000 to 50 copies per PCR of 
the P-globin gene using human placenta DNA (Sigma, art. no 
D701 1). The standard curve had a slope of —3.4, y intercept of 
39.9 and r 2 of 0.99. The PCR efficiency calculated from slope was 
97.5%. For the calculations of number of human cell per sample, 
the copy number of P-globin was divided by two. We assumed that 
each human cell carries two P-globin gene copies and that the 
diploid genome equivalent contains ~6.6 pg DNA. No-template 
controls of water samples were tested of both the HPV154 PCR 
and the human P-globin gene PCR. All samples were analyzed 
triplicate. Coefficient of variation was calculated for each triplicate 
measurement of viral copy number per human cell. In the 
quantitative PCR, negative controls showed no Ct values. 

Bioinformatic Analysis 

General sequence handling and feature identification were 
carried out using the UGENE software vl.11.5 (http: //genome. 
unipro.ru). Putative proteins from ORFs were generated by 
UGENE and they were then searched for similarities with other 
proteins using BLASTp (http://blast.ncbi.nlm.nih.gov). Proteins 
were analyzed for unique domains with ScanProsite (http:/ /www. 
expasy.ch/prosite) and SMART [47], including searches in the 
Pfam database (http://pfam.sanger.ac.se/search) [48]. A Python 
v3 script was made to compare pairwise identities from multiple 
sequence alignments (available upon request), in which all 
differences, including terminal gaps, were counted. 



Phylogenetic Analysis 

In order to avoid missing any gamma-related virus, all complete 
genomes were retrieved from GenBank (914 accessions), from 
which 723 unique genomes remained after dereplication using 
Usearch v6. 0.307 [49]. A workflow made using UGENE software 
[50] was used to extract the FAP region from LI ORFs, align it 
using uMUSCLE [50,51], and the resulting alignment was then 
converted to phylip format using ALTER website [52] . It was used 
later to make an ML-tree with RaxML (400 rapid bootstrap 
inferences and thereafter a thorough ML search). The branch 
containing all gamma-related viruses was selected using Dendro- 
scope v3.1.0 [53] (excluding beta types). From those sequences, LI 
was extracted, aligned, and dereplicated again to remove identical 
accessions or genomic variant sequences. The complete genomes 
of the remaining 91 unique sequences were manually edited to 
have LI at the 3' end and then they were aligned with uMUSCLE. 
It was later used to make a Maximum Likelihood phylogenetic 
inference using RAxML v 7.4.2 [54], compiled under Linux as 
AVX and PTHREADS versions, and run using the raxmlGUI 
vl.3. Gaps were treated as missing data and the General Time 
Reversible (GTR) under the gamma model of rate heterogeneity 
was selected as a nucleotide substitution model to make 20 
inferences with 150 thorough bootstrap replicates (automatically 
determined by the program using the majority rule tree based 
criteria, command "-N autoMR"). The sequence of the betapa- 
pillomaviruses HPV5 (acc. no. Ml 7463) and HPV9 (acc. no. 
X74464) was used to root the tree. A graphical representation of 
the tree was made with FigTree software v 1.3.1 [55] and edited 
with Inkscape v0.48. The GenBank accession numbers of the 
sequences in the tree are as follows: AsPVl, HQ625440; BPV11, 
AB543507; BPV12, JF834523; BPV3, AF486184; BPV4, X05817; 
BPV5, EU360723; BPV6, AJ620208; BPV7, DQ2 17793; BPV9, 
AB331650; CG2, JF966378; CG3, JF966379; ChPVl, 
DQ091200; CPV13, JX141478; CPV2, AY722648; CPV7, 
FJ492742; DdPVl, GU117620; FA69, KC108722; FD1, 
JF966375; FD2, JF966376; Fi864, KC311731; FS1, JF966373; 
HPV4, X70827; HPV48, U31789; HPV50, U31790; HPV60, 
U31792; HPV65, X70829; HPV88, EF467176; HPV95, 
AJ620210; HPV101, DQ080081; HPV103, DQ080078; 
HPV108, FM212639; HPV109, EU541441; HPV112, 
EU541442; HPV116, FJ804072; HPV119, GQ845441; 
HPV121, GQ845443; HPV123, GQ845445; HPV126, 
AB646346; HPV127, HM011570; HPV128, GU225708; 
HPV129, GU233853; HPV130, GU117630; HPV131, 
GUI 17631; HPV132, GUI 17632; HPV133, GUI 17633; 
HPV134, GUI 17634; HPV135, HM999987; HPV136, 
HM999988; HPV137, HM999989; HPV138, HM999990; 
HPV139, HM999991; HPV140, HM999992; HPV141, 
HM999993; HPV142, HM999994; HPV144, HM999996; 
HPV146, HM999998; HPV147, HM999999; HPV148, 
GU129016; HPV149, GUI 17629; HPV153, JN171845; 
HPV154, JN211193; HPV155, JF906559; HPV156, JX429973; 
HPV161, JX413109; HPV162, JX413108; HPV163, JX413107; 
HPV164, JX413106; HPV165, JX444072; HPV166, JX413104; 
HPV169, JX413105; HPV170, JX413110; KC5, JX444073; 
KN1, JF966371; KN2, JF966372; KN3, JF966374; McPV2, 
DQ664501; MmiPVl, DQ269468; MusPVl, GU808564; 
PphPVl, GUI 17621; PphPV2, GUI 17622; PsPVl, AJ238373; 
RnPVl, GQJ 80114; SD2, KC113191; SE87, KC108721; TtPVl, 
EU240894; TtPV2, AY956402; TtPV3, EU240895; TtPV4, 
JN709469; TtPV5, JN709470; TtPV6, JN709471; TtPV7, 
JN709472. 
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